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ABSTRACT 



We estimate flux ratios of the extrasolar planet X0-2b to its host star 
XO-2 at 3.6, 4.5, 5.8 and 8.0 micron with IRAC on the Spitzer Space Tele- 
scope to be 0.00081 ± 0.00017, 0.00098 ± 0.00020, 0.00167 ± 0.00036 and 
0.00133 ± 0.00049, respectively. The fluxes provide tentative evidence for a 
weak temperature inversion in the upper atmosphere, the precise nature of which 
would need to be confirmed by longer wavelength observations. X0-2b substellar 
flux of 0.76 X 10^ ergs cm~^ s~^ lies in the predicted transition region between 
atmospheres with and without upper atmospheric temperature inversion. 

Subject headings: stars:individual(X0-2) — binaries:eclipsing — infrared: stars 
— planetary systems 



1. Introduction 



The field of comparative planetology has burgeoned in the past year (jDemingj |2009| ): 



four exo-planet secondary eclipse observations have been published with all four IRAC 
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channels: HD 209458b ( 


^nutson et al. 


2008 


); HD 189733b ( 


Charbonneau et al. 


2008): XO- 


lb( 


Machalek et al. 


2008 


) and TrES-4 ( 


Knutson et al. 


2009). which were used to deduce the 



temperatures structure in the planetary upper atmospheres. Numerous scientific firsts were 
also announced: a r nodel dependent deriv ation of the radiative time constant of an extrasolar 
planet atmosphere (ILaughlin et al.l 120091 ) , first unambiguous d etection of water vapor in a 
Hot Jupiter atmosphere in HD 189733b (iGrillmair et al.l 120081). together with p rospects for 
detailed study of extrasolar planet atmospheres with JWST (ISeager et al.l 120081 ). 



Planetary stratospher es, which are common in the giant planets of the Solar System 
( jde Pater fc Lissauerll200ll ). form when high stellar fiux penetrates deep into the atmosphere 
where the atmosphere does not e fficiently radiate the heat away. Hot-Jupite r s were also ob 



serve d to possess stratospheres ( IKnutson et al.l l2008l : iMachalek et al.l l2008l : iKnutson et al. 



20091 ) and enhanced opacity at high altitude in the form of an extra absorber of optical 
l ight has been s ugges te d as the physic a l caus e of upper atmo s pheric te mperature inv e rsions 
faubenv et aP J2003h : burrows et aP J2008h : IPortnev et all hoO^ ). ISpiegel et all J2009h 
have recently argued against vanadium oxide (VO) as the extra optical absorber and sug- 
gested that even the previously favored gaseous optical absorber titanium oxide (TiO) would 
probably rain out of the upper atmosphere unless unusually high levels of macroscopic mix- 
ing exist (Eddy diffusion coefficients K^z ~ 10^-10 -^^ depending on particula te size of the 
TiO condensates) com pared to Jupiter ( K^ ^ ~ 10^ (jde Pater fc Lissaueiil200ll )) to overcome 
gravitational settling. IZahnle et al.l (|2009l ) has considered the role by which absorption of 
UV and visible light by S2 and S3 can lead to formation of thermal inversion in upper atmo- 
spheres of Hot-Jupiters and suggested planetary metallicity, in addition to substellar flux, 
as the determining factor for the presence of a temperature inversion. 



Hot- Jupiter atmosphere models (IBurrows et al.ll2007al . l2008l : ISpiegel et al.ll2009l : iFortney et al 



20081 ) show that water and CO opacity deflne the r=2/3 decoupling layers and the relative 
temperatures at those layers determine the relative brightness at 3.6 micron and 4.5 micron. 
In the presence of an optical absorber in the upper atmosphere the 4.5 micron planet flux 
is higher than the 3.6 micron flux and vice-versa for models without an optical absorber 
in the upper atmosphere, which serves as an important diagnostic for detection of thermal 
inversions in upper atmospheres of Hot-Jupiters (see Fig. [1] for the X0-2b Temperature / 
Pressure proflle). 

A detailed study of the IR secondary eclipse planetary spectra of HD 209458b, HD 



89733b, TrES-1, HD 149026b and non-eclipsing HD 179949b, and v And b bv Burrows et al. 



(I2OO8I ) suggests that the presence of such an upper atmospheric absorber might be dependent 



- 3 - 




Fig. 1. — (top) Te mperature / P r essure profile s for the atmo s phere of X0-2b following 
the methodology of iBurrows et al.l (l2007al l2008l ): ISpiegel et al.l (120091 ) for heat redistribu- 
tion parameter P„=0.3 with no upper atmospheric optical absorber (dot-dashed line) and 
a corresponding model with a uniform upper atmospheric absorber (solid line) (absorption 
coefficient Ke=0.1 cm^/g) with depths corresponding to emission in the IRAC channels de- 
noted. 



on the flux from the star at the sub-stellar point on the planet as well a s seco nd order effects 
like metallicity and planetary surface gravity. In the iBurrows et al.l (120081 ) interpretation 
planets with high sub-stellar point flux (e.g., HD 209458b, OGLE-Tr-56b, OGLE-Tr-132b, 
TrES-2b and X0-3b) would have extra optical absorber in the upper atmospheric layer and 
water features in emission while planets with lower fluxes (XO-lb, TrES-1, X0-2b and HD 
189733b) wou l d have no such extra absorber and would possess water features in absorption. 
Fortney et al.l (120081 ) also suggest a similar division of planetary spectra based on incident 
stellar flux. Observations of X0-2b and other planetary systems directly constrain and test 
the incident flux threshold necessary for the occurrence for a thermal inversion in the upper 
atmosphere of a Hot- Jupiter. 



Based on the p laneta ry sub-stellar point flux from the star, both IBurrows et al.l (120081 ) 
and iFortney et al.l (120081 ) predicted that XO-lb should not have exhibite d a thermal in 



version in its upper atmosphere, yet a temperature inversion is observed (IMachalek et al. 



20081 ). The incident flux on X0-2b from its parent star (0.76x 10^ ergs cm~^ s~^) lies in the 
predicted tra nsition region between upper atmospheres wi th or without thermal in versions as 
delineated by IBurrows et al.l (120081 ) ; iFortney et al.l (120081 ) . IBurrows et al.l (120081 ) predicted 
a "weak stratosphere" for X0-2b. 
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XO-2 (2MASS J07480647+5013328) has high metalhcity, [Fe/H]=0.45 ± 0.02, high 
proper mot ion, iJ,tn^=157 mas yr~^, and a common proper motion stellar companion with 31" 



separation JSurke et al.l boOTh. The planet X0-2b has an orbital period of 2.6158640 ± 0.0000016 days 
(IFernandez et al.ll2009l ) and orbits around the northern declination component of the XO-2 
stellar binary system (see Fig. [2]) . We presen t observations of the infrared spectral en- 
ergy distribution (SED) of the planet X0-2b (IBurke et allboOTI ) in all 4 IRAC channels 
obtained during secondary eclipses observed with Spitzer Space Telescope. By comparing 
X0-2b's SED with atmospheric models, we test for the presence of a thermal inversion in 
the upper atmosphere of X0-2b. 



Observations 



The InfraRed Array Camera (IRAC; iFazio et al.ll2004l ) has a field of view of 5.2' x 5.2' 
in each of its four bands. Two adjacent fields are imaged in pairs (3.6 and 5.8 microns; 4.5 
and 8.0 microns). The detector arrays each measure 256 x 256 pixels, with a pixel size of 
approximately 1.22" x 1.22". We have observed the X0-2b system in all 4 channels in 
two separate Astronomical Observing Requests (AORs) in two different sessions: the 3.6 
and 5.8 micron channels for 5.94 hours on UT 2007 November 16 (AOR 24462080) and the 
4.5 and 8.0 micron channels for 5.94 hours on UT 2007 November 19 (AOR 24462336). We 
used the full array 2s+2s/12s frame time in the stellar mode in which the 3.6 micron and 
4.5 micron bands are exposed for two consecutive 2s exposures while the 5.8 micron and 
8.0 micron bands are integrating for 12s to prevent detector saturation. Figure [2] shows a 
representative IRAC 3.6 micron image. 

We used the standard IRAC Basic Calibrated Data (BCD) products (version 16.1) 
described in the Spitzer Data Handbool!0, which includes dark frame subtraction, multiplexer 
bleed correction, detector linearization, and flat-fielding of the images. The starting point 
for our analysis were the BCD images. We converted the times recorded by the spacecraft 
in the FITS file header keyword DATE-OBS to heliocentric Julian dates using the orbital 
ephemeris of the spacecraft provided by the Horizons Ephemeris Syste ng. We flagged cosmic 
ray pixels in the images, which resulted in [0.8; 0.4; 2.4; 3.1] % of photometric points in the 
[3.6; 4.5; 5.8 and 8.0] micron time series to have a flagged pixel in the photometric aperture 
and which were not used in the analysis. We did not resample the pixels in any way during 
our analysis doing so could compromise the photometry. 



^http : //ssc . spitzer . caltech. edu/ irac/dh/ 



' http : //ssd . j pi . nasa . gov/ 
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Fig. 2. — The field-of-view of the XO-2 binary system in the IRAC 3.6 micron channel. 
North is up, East is to the left. X0-2N hosts the transiting extrasolar planet X0-2b. 



To determine the scalar background for the photometry we subtracted the zodiacal 
background in each channel by constructing a histogram of all pixels in each image and 
fitting a Gaussian to the distribution of the zodiacal background brightness. We evaluated 
the centroids of both North and South components of X0-2b with the IDL CNTRD routine 
which locates where the second derivatives of x and y pixel position reach 0. We have found 
that CNTRD centroids have lower rms in both x and y pixel directions than the centroids 
produced by IDL routine GCNTRD, which fits Gaussians to marginal x and y distributions. 
The pointing in the 3.6 micron channel varied by ~ 0.10 pixels in x and ~ 0.06 pixels peak- 
to-peak in y and by ~ 0.15 pixel in x and ~ 0.20 pixel in y in the 4.5 micron channel. The 
shifting of the stellar centroid within a pixel, which have sub-pixel sensitivity variations, 
resulted in a modulation of the stellar flux in the 3.6 micron and 4.5 micron channels (de- 
scribed below). 



2.1. 3.6 and 4.5 micron time series 

We performed aperture photometry on the 3.6 micron and 4.5 micron background sub- 
tracted images with aperture radii ranging from 2.5 to 5.0 pixels in 0.5 pixel increments for 
both components of the XO-2 binary. The sky level was elevated during the first hour of 
observations in the 3.6 micron channel (see Fig. [3]). This behavior is oppos ite the asymp 



totic behavior in the 3.6 micron channel for the similarly bright star TrES-4 (IKnutson et al. 



20091 ). The 3.6 micron and 4.5 micron time series exhibited a sharp increase during the 
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Fig. 3. — The average sky count rates in the annulus with radii between 13 and 50 pixels 
around the star X0-2N for the 3.6 micron IRAC channel observations. Notice the asymptotic 
increase of the sky count rates in the first ~ 1 hour of observation. Similar trend was observed 
for the other binary component X0-2S. 



first ~ 20 minutes of exposure for XO-2 and the 2 calibrators, presumably as a result of the 
instrument reaching a new equilibrium after previous observations. Such relaxation effects 
can reach several percent and usually stabilize within the first hour of observations of a new 
target. We have rejected the first 180 points in the 3.6 micron and 4.5 micron time series. 

The aperture radius was selected by choosing the lowest rms for out-of-eclipse points. 
Aperture of 2.5 pixels was chosen for the 3.6 micron time series, with an out-of-eclipse 
rms of 0.0044 after systematic effect removal which is 1.02 times the photon noise limit 
based on stellar brightness, background flux and detector read noise. Similarly, for the 
4.5 micron series an aperture of 2.5 pixels resulted in the lowest out-of-eclipse rms of 0.0063 
which 1.05 times the photon noise limit. We applied the appropriate aperture correction 
for each channel to the stellar flux value according to the Spitzer Data Handbook of [1.112, 
1.113, 1.125, 1.218] for the [3.6, 4.5, 5.8 and 8.0] micron channels, respectively. 

A strong correlation between the sub-pixel centroid and stellar brightness was observed 
in both the 3.6 micron and 4.5 micron c hannels, with flux amplitudes of ~ 1.3 % and .5 



2006; Machalek et al. 2008; Knutson et al. 



, respectively. This well s tudied effect (ICharbonneau et al.l l2005l : [Morales- Calderon et al. 



20091 ) is due to the InSb detector intrapixel sen- 
sitivity variations as the spacecraft jitters ~0.1 - 0.3 arcsec in orientation over a period of 
~3000 second^. The uncorrected sub-pixel intensity variations are clearly visible in the time 



|http : //ssc ■ spitzer . caltech.edu/docijments/exopleLnetmeino.txt 
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series of XO-2 in the 3.6 micron and 4.5 micron channels in Fig. HI 

Despite the similarity between the XO-2 stellar binary components (same brightness J 
= 9.74± 0.02 and similar color {J-il)NORTH = 0.40± 0.03 vs. iJ-li)souTH = 0.37± 0.03) the 
amplitude of their the sub-pixel sensitivity variations in the 3.6 micron channel for the North 
and South component of X0-2b is markedly different (1.3 % for X0-2N vs 0.7% for X0-2S, 
see Fig. H]). This is most likely due to the different positions of the two XO-2 components 
with resp ect to the edge o f the p ixel and hence being affected differently by the pixel response 
function. iLaughlin et al.l (120091 ) have noted in their 30-hour long 8.0 micron time series of 
the HD80606 and HD80607 binary components, which also have same brightness and similar 
color, that due to different sub-pixel positions the detector systematic ramps are different. 

We have corrected for the sub-pixel intensity variations by fitting a function of 4 variables 
to the time series of X0-2b : 



Isubpixei = bi+ b2X + h^y + bit, 



T. 



where bj are coefficients; x and y are subpixel centroids of the stellar flux and t is the 
time from the start of observation s in days. We hav e tried quadratic and cross terms for 
the subpixel centroids correction (jPesert et al.l l2009l ) but that has not improved our 
and hence were not used. We fit the secondary eclipse light curves using the formalism of 



Mandel &: Agoll (120021 ) with no stellar limb darke ning as appropriate for a secondary eclipse. 



We adopt stellar and planetary parameters from iBurke et al 



98+°-°3 



R, 



J2007h R 



R^ = 0.97 



+0.02 
-0.02 -^^0 5 



R(71 



88.9 



+0.7 



ephemeris from [Fernandez et al.l (120091 ): 



; degrees, and a = 0.0369 ± 0.0002 AU with updated 



TciE) = 2,454,466.884670(//JD) + £;(2.6158640 days) . 



(2) 



We have fitted the 4 baseline detrending parameters (a constant, linear x-position, linear 
y-position and linear time terms) concurrently with the depth of the eclipse AF and the 
timing of the centroid AT (total of 6 fitting parameters) using a Monte Carlo Markov Chain 
(MCMC) wit h 10^ ite r ations such that the ratio of jumps for each parameter was between 



20-40 % (see ICregorvl J2005h : iMarkwardtl toO^ ). The initial 20% iterations were rejected 



We have redone the entire analysis using the updated planetary and stellar parameters of lFernandez et al 



([20091) and our results do not change withing uncertainties. 
51 Rj„p = 71,492 km. 
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to remove arbitrary starting conditions and the final parameters were chosen as the median 
value from the posterior probability distribution for each parameter. The adopted eclipse 
depth AF and centroid timing AT in minutes from the expected secondary eclipse mid-center 
time for an assumed eccentricity of zero are reported in Table [T] with errors obtained from 
the symmetric 66.8 % contours around the median of the posterior probability distribution 
from the MCMC runs. Best-fit eclipse curves binned in 6.0-minute intervals are plotted in 
Fig. [5] and the eclipse parameters are listed in Table [H They are the channel wavelength, 
eclipse depth AF, eclipse mid-center time in HJD and the timing offset AT. 

We find that the 3.6 micron time series of X0-2b exhibits a linear flux trend with 
a slope of hi = -0.011% ± 0.005% per hour. This trend is consistently removed from our 
time series by in cluding a linea r time term b4 as one of the MCMC 6 parameter fit, as 



described above. iKnutson et al.l (120091 ) noted a linear trend in their 3.6 micron time series 
of TrES-4 with a slope of +0.030 ± 0.004 % per hour and attributed such linear trend to a 
previously unknown instrumental effect of the detector. The 4.5 micron time series of XO- 
2b exhibits a similar linear trend to the 3.6 micron time series with a slope of -0.010% ± 
0.004 % per hour, which is lower than a corresponding slope for TrES-4 of -0.020±0.003% per 



hour (IKnutson et al.ll2009l ) even though X0-2N is 11 times brighter. 



To assess the amount of non-Gaussian correlated red noise left in the time series after 
systematic effect removal we have performed the complete MCMC analysis on the other 
binary component X0-2S in the XO-2 binary system, which does not have a short-period 
Jovian mass orbiting it (which would have easily been detected by radial velocity observa- 
tions) and thus does not exhibit a secondary eclipse. The best fit eclipse curves of X0-2S are 
plotted on the right in Fig. [51 The 3.6 micron, 5.8 micron and 8.0 micron channels exhibit 
eclipse depths consistent with zero while the 4.5 micron channel shows an eclipse with depth 
0.00064 ± 0.00033 based on Gaussian statistics with a time shift of +11.1 ± 24.8 min. We 
exclude the possibility that the spurious eclipse at 4.5 micron around X0-2S could be due 
to a transiting super-Earth (corresponding to Rp=2.4R0 around X0-2S) by noting that the 
eclipse is not observed in the 8.0 micron time series of X0-2S, which allows for maximum 
depths of 0.00002 ± 0.00001 and since it was observed at the same time as the 4.5 micron time 
series we can reject the transiting super-Earth hypothesis at more than lO-a. 

The significant non-zero eclipse depth around the control star X0-2S which does not 
have short-period Jovian-mass planet suggests that the 4.5 micron time series for both X0-2N 
and X0-2S underestimates the amount of time correlated red-noise even after our calibration 
procedures. 
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2.2. 5.8 and 8.0 micron time series 

Aperture photometry was performed on the images with cosmic ray pixels flagged with 
an aperture radius between 2.5 and 6.0 pixels in 0.5 pixel increments. The optimal size of the 
aperture was determined by minimizing the rms scatter in the light curve for observations 
outside of the eclipse. Aperture of 3.5 pixels was chosen for the 5.8 micron time series, 
with an out-of-eclipse rms of 0.0066 which is 32 % greater than the the theoretical noise 
based on stellar brightness, background flux and detector read noise. Similarly for the 8.0 
micron series an aperture of 4.0 pixels resulted in the lowest out-of-eclipse rms of 0.0054 
which is 35% higher than the photon noise limit. The internal scattering of photons inside 
5.8 micron and 8.0 micron Si:As arrays is likely responsible for the fact that we do not 
approach the Poisson lir nit as closely as in the 3 .6 micron and 4.5 micron channelsfl Unlike for 
the XO-lb observations iMachalek et al.l (120081 ) . where the flrst ~ 30 minutes of observations 
were rejected as the instrument settled into a new equilibrium state, we do not observe such 
an initial ramp-up in the 5.8 micron and 8.0 micron time series and conversely we do not 
reject any data points from the beginning of the time series in these channels. 

Fig. m shows intensity variation with time, which is caused by changes i n the effective 



gain o f individual pixels over time. This effect has been observed before by iDeming et al. 



( I2OO5I ). both in the IRAC camera and in the IRS and MIPS 24 micron cameras and is depen- 
dent on the illumination level of the individual pixel (IKnutson et al.ll2007l . l2008l ). Pixels with 
high illumination will reach their equilibrium within ~1 hour, but lower illumination pixels 
increase in intensity over time, approximately proportional to the inverse of the logarithm 
of illumination. The detector ramp intensity in the 5.8 micron time series decreased in flux 
during the ~ 6 hours of observat ion by ~ 0.1 % for X 0-2N and by ~ 0.5 % for X0-2S, 



following a similar trend seen by iKnutson et al.l (120081 ) in their 5.8 micron time series of 
brighter HD 209458b. 

We detect a nonline ar increase in the bri ghtness of X0-2N and X0-2S in the 8.0 mi- 
cron channel similarly to IKnutson et al.l (120071 ) who have reported a nonlinear flux increase 
over time with the 8.0 micron IRAC detector. Most recently the gain variation with time 
and illumination has be en conflrmed in the extended duration (30 hours) observations of 
HD80606 and IID80607 (iLaughlin et al.ll2009l ). who found that even for similarly bright and 
similar color binary, the 5.8 micron and 8.0 micron time series "ramp" can be different, 
presumably due to different position of the stars in relation to the edge of the pixels. 

To properly account for the way in which the non-linear flux ramps affect our estimate 



http : //ssc . spitzer . caltech. edu/documents/iracjnemo . txt 
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of the secondary eclipse depth and timing we have performed MCMC simultaneous fitting 
of the 3 ramp parameters (see Eq. [3]) and the eclipse depth AF and centroid timing AT (5 
parameter fit overall): 

Imodei = ai + 02 X ln{At + 0.05) + as X {ln{At + 0.05))^ (3) 

where Imodei is the normalized model flux, At is the time in days since the beginning of obser- 
vations (the 0.05 factor is included to prevent singularity at At =0 ) and aj are coefficients. 

The MCMC fitting was performed in the same way as for the 3.6 micron and 4.5 mi- 
cron channels with uncertainties obtained from the symmetric 66.8 % contours around the 
median of the posterior probability distribution from the MCMC runs. The eclipse depth 
in the 5.8 micron channel of 0.00167 ± 0.00036 represents an overall 4.7-0" detection of the 
secondary eclipse based on the cumulative SNR of all the in-eclipse points and similarly the 
8.0 micron channel eclipse of 0.00133 ± 0.00049 has an overall SNR of 2.7. The resultant time 
series were normalized using out-of-eclipse points and binned into 6.0-minute bins (Fig. [5]) 
for viewing clarity. 



3. Discussion 

To check whether our results depend on the aperture radius we have redone the MCMC 
analysis for photometry with aperture radii between 2.5 and 4.5 pixels for the 3.6 micron and 
4.5 micron time series and aperture radii from 3.0 to 5.0 pixels in the 5.8 micron and 8.0 mi- 
cron time series and obtained consistent results for the eclipse depth AF and centroid timing 
AT within errors to our adopted values from Table [TJ 

To test the robustness of our data reduction and MCMC analysis technique and con- 
sistency with other observations in the IRA C full-array mod e we h ave reanalyzed the IRAC 



secondary eclipse time series of XO-1 from iMachalek et al.l (120081 ) in all 4 IRAC channels 



with the new MCMC pipeline. The eclipse depths AT of XO-lb with our MCMC pipeline 
are still in agreement with a model of a thermal inversion and an extra upper atmospheric 
absorber of uniform opacity of = 0.1 cm^/g and redistribution parameters of P„ = [0.3] at 
[1.8,0.8,0.4]-cr at 3.6 micron,4.5 micron and 5.8 micron bands, respectively. The secondary 
eclipse of XO-lb in the current MCMC analysis of 8.0 micron time series is 3.8-0" above 
the thermal inversion model above but clearly inconsistent at 16.1-0" with an atmospheric 



model without a temperature inversion. Our conclusion from iMachalek et al.l (120081 ) for a 
temperature inversion in the upper atmosphere of XO-lb caused by an optical absorber of 
uniform opacity of Kg = 0.1 cm^/g is thus reinforced. We have thus demonstrated that our 



- 11 - 



MCMC reduction and analysis pipeline is robust a nd the secondary eclips e depth estimates 
are consistent with previous pipeline versions from iMachalek et al.l (120081 ) . 



The eclipse mid-center timings for X0-2b in Table [T] are individually consistent within 
uncertainties wi t h zero timing residuals for a circular orbit based on the ephemeris by 
Fernandez et al.l (120091) and a c ombin ed mid-eclipse timing offset of -1.0 ±9.7 min. Us- 
ing the equation of (e.g. iKopall (119591 ) Eq. 9.23): 



e X cos^uj) 



nAt 



(4) 



where e is the eccentricity, uj is the longitude of periastron, P is the orbital period, and At is 
the centroid time shift from expected time of secondary eclipse, allows us to set a S-a upper 
limit on excos(u;) < 0.012. 

The X0-2b eclipse depths in Fig. |6] show several trends. The flux ratio of X0-2b to the 
star peaks in the 5.8 micron channel, with a decrease towards the 3.6 micron and 4.5 mi- 
cron channels and a slight decrease towards the 8.0 micron channel. The solid line and band 
averages represented as o pen squares in Fig. [6| depict a n atmospheric mode l of X0-2b, fol- 
lowing the methodology of iBurrows et al.l (l2007al . l2008l ): ISpiegel et al.l (l2009l ) with a thermal 
inversion and an extra upper atmospheric absorber of uniform opacity of Kg = 0.1 cm^/g 
and redistribution parameters of P„ = [0.3]. Pn=0 corresponds to no heat redistribution 
from the plane t ary d ay-side to the night-side and = 0.5 stands for full redistribution (see 
Burrows et al.l (120081 ) for details). The upper atmospheric temperature inversion model fits 
the observed eclipse depths within 1.0-cr in the 3.6 micron, 4.5 micron and 5.8 micron bands, 
respectively, while the 8.0 micron channel eclipse depth is offset by -lA-a from the model. 
The model without an upper atmospheric temperature inversion with the redistribution 
parameter P„ = [0.3] (dot-dashed curves) does not fit the observations by [ 1.6, 2.2]-cr, 
respectively, in the 4.5 micron and 5.8 micron channels while fitting the 3.6 micron and 
8.0 micron band eclipse depth within 1.0-cr. The eclipse depth fit in the 3.6 micron, 4.5 mi- 
cron and 5.8 micron bands thus provides evidence for a "weak" temperature inversion in the 
upper atmosphere of X0-2b. X0-2b would thus be the second Hot- Jupiter along with XO-lb 
with a sub-stellar point flux of less than ~1.0 x 10^ erg cm ~^ s to possess a temperature 
inversion in its upper atmosphere. 

The possibili t y of th e rmal inversion in a planetary upp e r atnaosphere has been suggested 
bv iHubenv et al.l J2003h : liro et all J2005h : burrows et aP J2006I . l2007ah and IPortnev et al 
( 120081 ). A thermal inversion in the planetary upper atmosphere ha s been invoked for in - 
terpretation of broadband IR spectra in the case of HD209458b (IKnutson et al.l |2008|); 
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XO-lb flMachalek et al.l 120081): HD149026b flHarrington et all 120071 ) and recently TrES-4 
( iKnutson et al.ll2009l ). The presence of an extra optical absorber in the upper atmosphere, 
of yet unknown composition, would yield a thermal inversion in the planetary upper atmo- 
sphere and the presence of the water features in emission. 



Burrows et al.l (120081 ) and iFortney et al.l (120081 ) suggested that the presence of the extra 



optical absorber in the upper atmosphere might be correlated with the incident flux from 
the star at the sub-stellar point on the planet. X0-2b substellar flux of (0.76 x 10^ ergs 
cm~^ s~^) lies in the transition region between atmospheres with or without upper atmo- 
spheric temperature inversions. XO-lb which has an even lower substellar flux of 0.49 x 
10^ ergs cm~^ s~^ yet possesses unambiguous signs of a temperatur e inversion in the u pper 



atmosphere, which suggests that secondary effects like metallicity (jZahnle et al.ll2009l ) and 
planetary surface gravity may also determine the nature of the temperature inversion of the 
upper atmospheres of Hot-Jupiters. 



Confirmation of the weak temperature inversion of X0-2b would be possible with longer 
wavelength observations because the contrast between models with or without a thermal in- 
version is high (see Fig. [6]). Further study of Hot- Jupiter upper atmospheres, especially of 
planets that lie in the predicted transition region of the substellar flux (~ 0.6 - 1.0 x 10^ ergs 
cm-2 s"^) like HAT-P-1, WASP-2b, HD197286, should reflne the sub-stellar flux boundary 
with respect to the presence/absence of an upper atmospheric thermal inversion. 



4. Conclusion 

We report the estimated flux ratios of the planet X0-2b in the Spitzer Space Tele- 
scope IRAC 3.6, 4.5, 5.8 and 8.0 micron channels. The fluxes are consistent with a weak 
temperature inversion in the upper atmosphere of X0-2b. The atmospheric model with 
an upper atmospheric temperature inversion with an extra optical absorber of opacity of 
Ke = 0.1 cm^/g and redistribution parameter P„=0.3 provides the best flt to the data in 
the 3.6 micron, 4.5 micron and 5.8 micron channels with a mild -lA-a inconsistency in the 
8.0 micron channel. 

The presence or absence of the stratospheric absorber and thermal inversion layer has 
been linked to the flux from the parent star at the sub-stellar point on the planet. The XO- 
2b sub-stellar point flux of ~0.76 x 10^ erg cm ~^ s ~^ is the second lowest so far reported 
for a planetary atmosphere with a thermal inversion. 
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Table 1. Secondary eclipse best fit parameters 



A 


Eclipse Depth AF 


Eclipse Center Time 


Time offset AT 


(microns) 




(HJD) 


(min) 


3.6 


0.00081 ± 0.00017 


2454421.10118 ± 0.01675 


-6.9 ± 24.1 


4.5 


0.00098 ± 0.00020 


2454423.72707 ± 0.01306 


7.8 ± 18.8 


5.8 


0.00167 ± 0.00036 


2454421.10719 ± 0.01315 


1.8 ± 18.9 


8.0 


0.00133 ± 0.00049 


2454423.71817 ± 0.01221 


-5.1 ± 17.6 
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Fig. 4. — (left) Secondary eclipse observations of X0-2b with IRAC on Spitzer Space Telescope in 
3.6 micron, 4.5 micron, 5.8 micron and 8.0 micron channels (from top to bottom) binned in 6.0- 
minute intervals and normalized to 1 and offset for clarity. X0-2b orbits X0-2N the northern 
component of XO-2 binary, (right) The light curve of the southern component X0-2S of the binary 
system (with no planet) plotted to the same scale and included for comparison. The overplotted 
solid lines do not represent a fit to the time series, but rather show the corrections for the detector 
effects (see text). 
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Fig. 5. — (left) Secondary eclipse of X0-2b around the binary component X0-2N observed 
with IRAC on Spitzer Space Telescope in 3.6, 4.5, 5.8, and 8.0 micron channels (top to 
bottom) corrected for detector effects, normalized and binned in 6.0-minute intervals and 
offset for clarity. The best-fit eclipse curves are overplotted. (right) The time series of the 
southern component X0-2S of the XO-2 binary, which has no eclipsing planet, after detector 
correction. The non-zero depth in the 4.5 micron channel of X0-2S is a spurious residual of 
systematic effects that mimics an eclipse. See §2.11 for details. 
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Fig. 6. — Spitzer Space Telescope IRAC secondary eclipse depths for X0-2 b with MCMC 
error bars (filled squares) . The predicted emission spectrum of the planet (jSurrows et al. 



2nn7al . bnnsi : ISpiegel et aLlboOoh with an upper atmospheric absorber of = 0.1 cm^/g and 



redistribution parameter of Pn=[0.3] is plotted as a solid line. A model with no atmospheric 
absorber and a redistribution parameters of P„=[0.3] is over plotted with dot-dashed line 
(see [|3] for details). The band-averaged fiux ratios are plotted as open squares and open 
circles for the models with and without an extra upper atmospheric absorber, respectively. 
The normalized Spitzer Space Telescope IRAC response curves for the 3.6-, 4.5-, 5.8-, and 
8.0 micron channels are plotted at the bottom of the figure (dotted lines). 



